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In this paper certain type of biometric measurements has been used to identify 
the cone beam computed tomography (CBCT) radiograph of the subject in a 
fast and reliable way. Where the CBCT radiograph of a person is used as a 
data and stored in database for later use in a person’s recognition process. The 
aim of this research is to use various stages of the preprocessing operations of 


the CBCT radiograph to obtain the clearest possible image that will help us in 
the identification process more easily and precisely. The contourlet 
Keywords: transformation was used for feature extraction of each particular CBCT image 
CBCT and the results were processed by a new hybrid particle swarm optimization 
(PSO) named "contourlet PSO" algorithm (CPSO), which is faster and 
Contourlet produce more precise (due to apply contourlet algorithm) than traditional 
CPSO PSO. The proposed algorithm (CPSO) gave a detection ratio of 98% after its 
Directional filter application on 100 CBCT radiographs. 
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1. INTRODUCTION 

The personal secret identification number can be lost, forgotten or simply difficult to memorize. Also 
it may be stolen or hacked on in many occasions. Therefore, biometrics can be used for identity access as an 
alternative [1], [2]. The main advantage of biometric features is that these are not prone to theft and loss, and 
do not rely on the memory of their users. Moreover, biometrics, such as finger prints, iris print, hand geometry, 
ear shape, face and teeth scans don’t change significantly over time and it is a difficult for a person to alter his 
own physiological biometrics or imitate other individuals. So, the biometric record took a lot of interest in the 
last decades as a safer method for personal identification. Each one of these biometric records has its advantages 
and disadvantages [1]. Dental biometrics uses dental structure information for the automatic identification of 
human remains [3]. So the aim of this article is to present a medical system for personal identification 
depending on forensic dentistry where an artificial intelligent system have been made using the cone beam 
computed tomography (CBCT) for personal identification, and despite the fact that teeth shape can be subjected 
to change over time, but the shape, size of jaw and facial bones usually are still clear and constant without any 
change, so these features can be efficiently used for personal identification and in artificial intelligence 
algorithms. for this reason we chose CBCT images radiographs as a database for the proposed system. The rest 
of this paper is outlined as shown: section 2 deals with related works, section 3 deals with CBCT, section 4 
deals with contourlet algorithm, section 5 deals with the proposed algorithm, results and analysis are shown in 
section 6, section 7 deals with conclusion and future works, finally, the acknowledgement. 

Dighe and Shriram introduced a method for identifying the individual identity from dental 
information.The proposed method included three main steps: feature extraction using morphological operations 


Journal homepage: http://ijai.iaescore.com 


398 0 ISSN: 2252-8938 
(mathematical morphology), dental code generation, and dental code matching.The method was implemented 
on 30 X-ray database [4]. Karunya et al. [5] proposed a new system which consists of two stages: first stage 
includes extraction of features and the second stage is matching. In the first stage the contour model was used 
for contour extraction. The second stage includes two sub stages evaluation the distances of the images and 
identification of the subject. The system was applied on ten normal images and fifty-five orthopantomagram 
(OPG) images. Rehman et al. [6] presented an efficient method for human authentication correctly which 
consists of five main processing stages: preprocessing, segmentation, processing steps for segmentation, 
feature extraction and biometric analysis. The method was tested on colored teeth images for 14 persons and 
dental radiographs images for 45 persons. 

Oktay et al. [7] presented a method for distinguishing humans by comparing two-dimension 
panoramic dental X-ray images. First each tooth is detected and labeled using support vector machine and 
graphical probabilistic models. The matching ratings between images were calculated based on an appearance 
of the tooth and the geometric similarities. Shaker et al [8] introduced a method for identity identification based 
on X-ray image. The method included three stages: preprocessing, feature extraction, and matching. The 
preprocessing wasmean, distance and standard derivation (STD), feature extraction was variance and principal 
component analysis (PCA). The best results showed an identification rate of 89%. Khudhur et al. [9] introduced 
a system to construct a database containing dental ante mortem radiographic features, later used for post 
mortem dental matching. The algorithm which applied on X-Ray image included three stages: segmentation of 
images, classification and extracting features. These features were STD, euler number and area taken from 
bite-wing X-ray image. 

In this research a hybrid method was proposed for discrimination CBCT radiation images. This 
method is a combination of discriminative features of contourlet coefficients and intelligent features of particle 
swarm optimization (PSO) algorithm. Firstly, The CBCT radiographic images was preprocessed through 
different steps to obtain clearest possible image That make the identification process simpler and more reliable. 
Then contourlet transformation was used for feature extraction of each particular CBCT image. Finally, the 
PSO algorithm was implemented on the extracted features for identification process. The new hybrid PSO 
method was faster and yielded more accurate result, and the use of CBCT radiographs images which gives 
information not found with the traditional two-dimension imaging added a strength point to the research. Also, 
the process of hybridizing the pso algorithm with a contourlet transformation added strength to the algorithm 
which improved the efficiency of the algorithm. A comparison with related previous studies were described in 
Table 1. 


Table 1. Comparisons with related studies 


No Research name Algorithms/methods No. of images Recognition rate 
1 Dental biometrics for human mathematical morphological 30 X-ray images 90% 
identification based on dental operations 
work and image properties in 
Periapical radiographs/ 2012 
2 Human identification using Shape registration method, euclidian ten normal images 72% 
dental biometrics/2014 distance and fifty five OPG 
images 
3 Human identification using Different methods for segmentation colored teeth Equal error rate (EER) 


dental biometric analysis/2015 and feature extraction images for 14 88.8% for colored images 


persons and dental 85.7% for dental 
radiographs images __ radiographs. 
for 45 persons 
4 Dental X-Ray based human — Standard deviation (STD), Euler 80 I 80 X-Ray 10% 
identification system for number & area taken from bite-wing images 
forensic/ 2017 X-ray image. 


Human identification with dental 
panoramic radiographic images 
/2018 

based 


Identification dental 


image/ 2018 


The use of contourlet 
transformations in hybridization an 
development of an intelligent PSO 
algorithm to distinguish cbct 
radiation image/ 2020 


support vector machine & 
graphical probabilistic models 


Mean, distance, standard derivation 
(STD), variance & principal 
component analysis (PCA) 


Particle swarm optimization & 
contourlet transformation 


206 X-Ray images 
of 170 various 
subjects 


75 75 x-ray images 
belonging to 115 
ppersons (five for 
every person 

100 CBCT 
radiographs images 


rank-1 precision of 
81% rank-2 
accuracy of 89% 


89% 


98% 
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2. THE CONE BEAM COMPUTED TOMOGRAPHY 

Arai et al. [10] and Mozzo et al. [11] working separately, presented the CBCT for the oral and 
maxillofacial applications and like computed tomography (CT), offered 3D investigation and increasingly 
precise imaging contrasted with 2D imaging. The financially savvy innovation of CBCT, prompted quick 
entrance into the field of dentistry with interest for responsibility of dental experts and dental instructors to 
investigate the uses of CBCT innovation. Radiographic assessment is necessary in diagnosis treatment planning 
in dentistry. Aside from packing three-dimensional life structures of the zone being radiographed into a two- 
dimensional picture, 2D imaging has many important drawbacks (including magnification, distortion, and 
superimposition), together prompting distortion of structures [12]. The applications of CBCT in dentistry 
include as: i) implantology: missing teeth substitution by dental implants requests precise visualization of the 
surgical site for the successful implant installation and to keep away from damage to adjacent important 
structures; ii) oral and maxillofacial surgery, orthodontics, end odontics, periodontics; iii) applications in 
temporom and ibular joint disorders; iv) applications in forensic dentistry: one of the parts of forensic dentistry 
is age estimation. Enamel is usually resistant to changes beyond ordinary wear and tear; on the other hand, the 
pulpodentinal complex displays physiologic and pathologicalalterations with aging. usually, to measure these 
changes, extraction and segmenting of teeth is vital, which isn't constantly a practicable decision. CBCT, 
conversely, provides a non-invasive substitute; and v) virtual treatment planning and simulations [13]. 


3. CONTOURLET ALGORITHM 

It is a true way to represent two-dimensional images, and is a new way to effectively represent the 
contour and texture of images [14]. The transformation consists of two-layer filters, where the laplace pyramid 
transformation is used to achieve multi-domain analysis and obtain discontinuous points. After this, the 
multidirectional analysis is carried out by the directional filter bank in order to connect the non-continuous 
points in the form of a linear structure [15], [16]. By incorporating the laplace pyramid and the directional filter 
bank it produces a multi-directional filter [17]. 


3.1. Laplace's pyramid 

Offers the means to achieve multiscale decomposition. In each step of decomposition it produces a 
lowpass downsampled version of the original image and a bandpass image. A coarse image with low 
frequencies and a more accurate image with additional high frequencies including point discontinuities are 
produced. This pattern can be repeated continuously in the lowpass image and is restricted only from the size 
of the original image Because of the downsampling [14]. 


3.2. Directional filter bank 

The directional filter bank (DFB) is designed to obtain high frequency contents such as smooth 
circumference and directional edges [18]. The DFB analyzes each detailed sub-range from laplace's pyramid 
(LP) to a number of directional sub-ranges. The package passing images from the LP are fed into DFB so that 
directional information can be obtained. The scheme of the multilayer decomposition (Contourlet). The 
merging between LP and DFB forms a dual filter bank which is called a pyramidal directional filter bank that 
analyzes the image into directional subdomains with multiple scales [14], [19]. Figure 1 shows the contourlet 
transform diagram. 


CBCTImage --4Y -—  # & =" 


Band pass of 
directional 


sub bands 


Figure 1. Contourlet transform diagram 
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4. PROPOSED ALGORITHM 
4.1. Preprocessing 

This is very important step which includes the basic steps performed on CBCT images. These steps 
improve the image quality and facilitate to use in the next stages. The steps of primary processing include the 

following Figure 2: 

— Taking CBCT images for a group of persons and incorporating them into the suggested database system 
as shown Figure 2(a) Then all the CBCT images were read from the available system database so the 
preprocessing starts for every image according to the following steps. 

— The CBCT images are converted from (RGB) to gray scale as shown Figure 2(b). 

— Improving the gray scale image color distribution by changing color pixels values through the use of 
contrast limited at adaptive histogram equalization (CLAHS), which is one of the morphological 
algorithms where the result used as a mask for the next [20]. This step improves the image to a great 
extent so the teeth, facial bones and the bones surrounding the teeth appear in dark colors with a higher 
contrast as shown in Figure 2(c). 

— Using another morphological algorithm known as (Erode) and applying it to the previous mask in order 
to restore the possible lost parts of the image in the previous stages through the use of SE=1, and “disk” 
function where they are used in rapidly restoring the lost parts without affecting the basic features of the 
image [21]. so, the teeth, jaws and facial bones become more clear. This step will result to a dark image 
that focuses on the cavities and bones as shown in Figure 2(d). 

— Applying one of the morphological methods named (morphological reconstruction) between the mask 
and image resulted from the previous stage where the contacted points in the images are extracted and 
rearranged in a new image [22]. This method utilizes the gray colors in a hybrid way and rebuilds them 
in a better condition. The results of this stage can be noticed in Figure 2(e). 

— The result from this stage gave a clear, less noise image which concentrating on the teeth and bones. the 
morphological methods and filters are considered effective in the image improvement process and color 
redistribution [23]. 

— Implementing sharpening for image improvement in order to make the borders be marked clearly 
including the edges of the teeth, bones and cavities. Also, dental fillings and the fixed dentures are seen 
in a clear white color as shown in Figure 2(f). 

— The image is converted to binary as shown in the Figure 2(g). 

— The black rows and colors surrounding the teeths and bones are deleted from the CBCT images as shown 
in Figure 2(h). 


Figure 2. Preprocessing steps results for CBCT image, where (a) the source image of CBCT, (b) the CBCT 
image after converted from (RGB) to gray scale, (c) the CBCT image after using morphological algorithms 
(CLAHS), (d) the CBCT image after applying morphological algorithm (Erode), (e) the CBCT image after 
applying morphological reconstruction, (f) implementing sharpening for the CBCT image, 
(g) the previous CBCT image converted to binary, and (h) after deleting the black rows and colors 
surrounding the teeths and bones from output images 
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4.2. Feature extraction 

The enhanced images resulted from previous stage are fed into 2-level contourlet transformation for 
feature extraction. The coefficients were unique for each individual image which gave a distinct feature for the 
particular image and thereby helped in the next step for recognition. 


4.3. Apply contourlet particle swarm optimization algorithm 

The main idea is to imitate the behavior of animals looking for food, such as fish, birds, or bees [24]. 
A variety of simple variations have been developed in order to increase the speed of convergence and the 
consistency of the solution found by the PSO [25]. PSO is developed for distinguishing images entered easily 
using extracted contourlet coefficients as the first generation in PSO algorithm. 

The use of contourlet coefficients resulting from contourlet transformations as input to the PSO 
algorithm made the algorithm faster to reach the solution and gave more accurate results, because the inputs to 
the algorithm are coefficients of a size less than the size of the original images, and therefore the number of 
iteration needed to reach the required results became less by 1/10, Where the optimal solution is continuously 
searched for in these transactions until it reaches the ideal solution depending on search and repetition. The 
PSO algorithm considers each solution as a particle and has two significant characteristics: position n and 
velocity vi. The two characteristics are related to each particle. So that n= (njl, nj2..., njN) and vi= (vil, 
vij2,...,vijN), where N reflects the dimensions of the problem and at each stage, the particles in the swarm are 
given a fitness function at each stage of the search for a solution. The speed and location values are modified 
in accordance with the (1) and (2): 


v=w*v-+t+cl ml1(n Best — 1) + c2m2(p Best — 1) (1) 
n=n+vAt (2) 


where W: represents weight of inertia responsible for regulating the impact of particles on past velocities; c1, 
C2: positive constants that are referred to as parameters of acceleration; m1, m2: random values in each 
appearance take on new values; At: represents the time steps; Best: is the best current position that has entered 
the particle or passed it until the present moment; pBestIt: is the best current location reached or moved by a 
particle of neighboring particles until the present moment. The contourlet PSO (CPSO) algorithm pseudo code 
will be: 


Input: Create community member sites by initialized position randomly of the practices: nj; 
(0) and velocityv,; (0). 

Output: - best position of the global optima n*. 

F(n;)=fitness, Which is calculated from equation [26]: 

SSIM (a, b) = 


((2u_a w_b + X_1)(20_ab + X_2))/((u_a*2 + p_b*2 + X_1)(p_a*2 + p_b*2 + X_2)) (3) 


Begin 

Repeat while max number of iteration is not reached do 
Begin 

For j=1 to number of particles 

IF F(n3) <F (nyena) then 

Nj = Nnew 57 

Update vj: using (1) 
Update nj: using (2) 
jt++; 

End 

End 


While fitness function for each member of the primary community is done using the similarity scale 
function (SSIM), where: pa: average of a, lb: average of b, u2 : variance of a, UZ : variance of b, Gap: 
Covariance of a and b, X;,X2,are Two variables to allow consistency in the partition process with a non-strong 
denominator. The following Figure 3 summaries the whole work. 


5. RESULTS AND ANALYSIS 

The results obtained after processing 30 different CBCT images of different people using the proposed 
algorithm (CPSO), as shown in Table 2 and Figure 4. It showed 100% detection rate in the training stage, it is 
the optimum value obtained from any recognition system. In the testing stage, another 57 different radiographs 
have been processed and the results were 98%. when minor changes to the radiographs were made on 13 images 
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belonging to the same persons these changes include (a change in the teeth such as extraction, the placement 
of fillings, or a change in the dimensions of the radiographs of the same person) The result was 100%. 


Figure 3. The work stages 


Table 2. Comparisons with related studies 


No. of processing _Type of processing Type of dataset No. of images DR ER WR 
1 Training familiar 30 100 0 0 
2 Testing familiar 57 98 1.8 5 
3 Testing Unfamiliar (after changes) 13 100 0 0 

120 
100 
80 @ No. of images 
” 
@ 
Fi 60 + DR 
40 AER 
20 > WR 
0 * & *% 
0 1 2 3 4 
No. of processing 


Figure 4. Results of proposed algorithm (CPSO) 


These results showed the strength of the system in identifying dental rays through the use of primary 
treatment, which showed the jaws very clearly, with the cancellation of all unwanted parts. Also, these high 
percentages show how important it is to use the contourlet method in the process of extracting features and 
identifying important points from the x-ray images and using them as inputs to the proposed algorithm and 
obtain matching results in properties by means of the SSIM equation, which was used in developing PSO 
algorithm and reaching the goals directly with the least mistakes. 

The scale used to calculate the overall degrees of discrimination in the system is detection rate (4). 
Where this measure represents the ability of the system to identify the person to whom the inserted dental rays 
belong. Therefore, it is calculated from the ratio of the number of correctly identified dental x-rays to the total 
number of x-ray images available in the system dataset. Therefore, the higher percentage of this scale is 
considered as a greater force for the system to discover the identity of the people to whom these rays belong. 
Thus, the value of detection rate is calculated as follows [26]. 


Detection Rate (DR) = (no. of correctly detected sampels) /(total no. of samples) * 100 (4) 


The extent of the system's errors in identifying people who have dental radiation inserted can be 
calculated by the scale Error Rate, and that was 1.8%, which is a small percentage that was calculated by 
the (5). This scale the lower its value, the higher the system has the ability to correctly achieve goals, due to 
the strength of the algorithm to find the right targets and avoid alien targets, This percentage can be calculated 
by the ratio of the number of images that were not found by the system from among the images entered to the 
number of images in the system dataset [26]. 


Error Rate = (no. of false detected samples) /(total no. of samples) * 100 (5) 
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Also, the wrong acceptance rate (6) in this system was 5%, meaning that any radiograph that does not exist in 
the radiograph dataset of the systemis not accepted. And this was done according to the following formula [27], 
[28]: 


Wrong acceptance rate = (no. of samples accepted error)/(total no.of images) *100 (6) 


6. CONCLUSION 

The initial treatment stages that were performed on the images contributed to the clear visibility of 
the radiograph and in turn helped to obtain high results in the identification process. using contourlet with the 
intelligent algorithm increases the ability of the algorithm to find the optimal solution effectively because it 
gives distinctive values for each image, which represent the characteristics and properties of that particular 
image. The challenges faced were the limited number of images available in the database which were used in 
the training and testing stages and difficulty distinguishing in the case of teeth falling out due to accidents or 
in the case of dental implant. The process of distinguishing dental rays is a good way to distinguish the identity 
of people, especially after physical changes have occurred to the external features of the human body over 
time. The radiograph identification rate with this system was high comparing with previous related studies, and 
thus the hybridization of the PSO algorithm with the Contourlet is considered a successful hybridization. 


7. FUTURE WORKS 

The possibility of making changes to a number of parameters of the proposed algorithm. This includes 
increasing the number of images acquired by the CBCT rays, also increasing the number of levels in the 
contourlet transformation, which leads to an increase in the number of features (coefficients) in the images and 
observe the magnitude of these effects on the results. Comparison of results using the proposed algorithm on 
both conventional and CBCT radiographs. Doing some noise on the images acquired in the database, applying 
to them the steps of the proposed algorithm, and observe the effect of this on the results. Use dental images to 
determine the identity of the deceased person by matching dental images of the deceased person provided by 
his relatives. The use of other methods in the discrimination process, such as machine learning. 
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